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Разработан аналитический метод г ‘идродинамического расчета радиального подшипника с адаптированным 
профилем его опорной поверхности. Даны оптимальные интервалы изменения безразмерных параметров, 
присущих сжимаемым микрополярным смазкам, а также параметра, характеризующего адаптированный 
профиль опорной поверхности подшипника. 
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Введение. Как известно, в качестве модели гидродинамической смазки в подшипниках скольже- 
ния в настоящее время широко используется сжимаемая микрополярная смазка. Естественно воз- 
никает необходимость не только аналитического прогнозирования оптимальных по несущей спо- 
собности и силе трения значений характеристик, присущих сжимаемым микрополярным жидко- 
стям, но и прогнозирование профиля опорной поверхности подшипника. 

Основной целью данной работы является разработка аналитического метода расчета радиаль- 
ного подшипника с адаптированным профилем его опорной поверхности, работающего на сжи- 
маемой микрополярной смазке. Оценка влияния параметра сжимаемости, параметров, присущих 
микрополярным смазкам, а также параметра, характеризующего адаптированный профиль опор- 
ной поверхности подшипника на основные рабочие характеристики подшипника. 

Постановка задачи. Рассматривается установившееся движение сжимаемой микрополярной 
смазки в зазоре радиального подшипника с адаптированным профилем его опорной поверхности. 
В полярной системе координат с полюсом в центре вала, уравнения контуров вала и кругового 
подшипника записываются в виде (рис.1) 


ЕК, "=и+ес0$0-азт®’0, (1) 
гДе о — радиус вала; ”, — радиус кругового вкладыша; е — эксцентриситет; а и © соответственно 
амплитуда и частота контурных возмущений. 





Рис.1. Схематическое изображение шипа в подшипнике с круговой опорной поверхностью: 
Со — контур вала; с! — круговой контур вкладыша 
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Более точное изображение адаптированного некругового контура опорной поверхности 
можно привести после определения оптимального (по несущей способности) значению параметра 
«. В дальнейшем предполагается, что е и а одного порядка малости. За исходные параметры бе- 
рется система безразмерных уравнений движения сжимаемой микрополярной смазки для «тонко- 
го слоя», уравнение неразрывности и уравнение состояния [1] 

0би оду Тр Эду у 10и 
ть — Е : =0 
д’ 07 ^а0 0’ М Ма 
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Эри дро 
оч =0, р=Ер. 2 
а о 
Размерные величины и’, о’, г', р’, у’, р’ связаны с безразмерными п, ъ, г, р, у, р соотноше- 
НИЯмМи 
о 
и'=Оки, о'=600, р'=р.р, у'= Ау, 
25 
р'=рр, р. еее. ‚п=и+аг, бЕИ-и. 
лок, 
Здесь 
2 2 
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А= и. 2. х , а и. ВА (3) 
р.:9 2+5 5 К 4 


где Л — параметр сжимаемости; и’о’ -— компоненты вектора скорости; р’ — гидродинамическое 
давление в смазочном слое; у — вектор скорости микровращения; и — динамический коэф- 
фициент вязкости для ньютоновской сжимаемой смазки; 7, к — коэффициенты вязкости 
сжимаемой микрополярной смазки; р’ — плотность; © — угловая скорость вращения вала; 
р’ — плотность смазки; р» -давление подачи смазки в случае подшипника конечной длины. 

Поскольку скорость движения направляющей считается достаточно большой, а поверхно- 
сти рассматриваемой пары трения являются шероховатыми, то ситуация соответствует так назы- 
ваемой квадратичной области течения жидкости, в которой потери давления на трение пропор- 
ционально квадрату скорости (формула Вейсбаха-Дарси) [2] 
ло? 

2 
где ^, — коэффициент потерь на трение, определяется экспериментально. 

Как видно из (2) и (3), помимо обычных безразмерных параметров, встречающихся в тео- 
рии сжимаемой смазки, вводятся и другие параметры для сжимаемой микрополярной смазки. Па- 
раметр взаимодействия № стремится к нулю при х>>0 и система уравнений (2) сводится к обыч- 
ным уравнениям Навье-Стокса для сжимаемой смазки. Параметр / имеет размерность длины и его 
можно трактовать как некоторую характеристику, зависящую от размера молекул смазки. Пола- 
гают, что микрополярные эффекты должны играть значительную роль, чем больше значение [, по 
сравнению с минимальной толщиной слоя смазки. Если последняя считается фиксированной, то 
чем больше / (т.е. чем меньше №'), тем ярче будут выражены микрополярные эффекты. 


Система уравнений (2) решается при следующих граничных условиях 


! 


В = 


! 


р 





и=Ьо=0, ео при х=0; р(0) = р(2п) =Ъ 


п 
и=0, о =0, у=0 при г=й(х) =1+1с0$0 — п, т 00, (4) 
где п=е/д, п =а/д, р, - давление питания в случае подшипника конечной длины. 
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Так как для сжимаемой микрополярной смазки 1/№,<<1, то пренебрегая в системе (2) чле- 
нами порядка 0(1/№\), будем иметь 











2 и 
и: Т@ ду —0, РР В (5) 
ди д" Л 40’ д? ду 09 
Интегрируя второе уравнение системы (5), с учетом граничных условий (4), будем иметь 
оо (6) 
170) п, 


Интегрируя первое уравнение системы (5), с учетом (4) и (6), получим 
И 1 ар! №5 , Гарй" г 
ъй л402 к А а02 й 
В случае уравнений (5), интегрируя по г уравнение неразрывности от 0 до й(0) с учетом 
(7), приходим к следующему уравнению 
3 
а р” рар и ый 
00| 12^ 40 6 к 


Интегрируя это уравнение, будем иметь 
Яр 12^ р 1 б 
аа, № |6 | (8) 
90 йЙр 26 к 
При больших угловых скоростях вращения вала, таких как О-—>(Л-—›), градиент гидро- 
динамического давления остается ограниченным при условии 


р ба 
№ —й > 0. 9 
2 . }-< (9) 


Отсюда, в рассматриваемом экстремальном случае для гидродинамического давления 





(7) 











б 
(с точностью до членов ов з 8, ), будем иметь 
И 


2 1 б 
ВЕНЫ 1- —(1с0$0-—п, $100) |, А — (10) 
28-п 28 в. 
Рассмотрим случай промежуточных значений угловой скорости вращения вала (фо, 


А-—>оо). 
Точное автомодельное решение системы (5), удовлетворяющее граничным условиям (4), 
будем искать в виде [3] 


ри 59% 00), ро О Г(",6), У =У(8), 
д+ 09 


у=(8), И=и(Е,0), о 








РР ео а(8) с —, &= с, = с0п8. (11) 
А 460 рк Й р т. 
Подставляя (11) в (5) и (4), получим 
3. — 
о (12) 
45 Е? 0 бх °&Е 


У =0 при &=0,Е=Е и(Е,0)=Р, 5=0 при Е=0; 
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и(Е,0)=0, 5=0 при &=5; [и(Е,0)4Е=0, 


0 


Зоо и 26.0 АЕ (13) 
Ра п 
Интегрируя систему (12) с учетом граничных условий (13), будем иметь 
м 6, =. 26 а -Е? [с 
у == (2-2), ®=—(-8), Кобе (14) 
7. 2 2 


где с =бр, а с, в дальнейшем определится из условия р(0)= р(2п) = Ро. Для определения гид- 


родинамического давления приходим к следующему нелинейному дифференциальному 
уравнению первого порядка 
а 5 ре 
Е 
д 49 п № 


Уравнение (14) решается численно при следующих граничных условиях р=2%, при 
Ра 


0=0, 9=2ж. Компоненты поддерживающей силы, сила трения, а также расход смазки опреде- 
ляются выражениями 


2п 2п 
В, =-пр, | рс0з046, К, =-пр, | рэ 646, 
0 0 





рр 
21° | 


В экстремальном случае (т.е. при Л->0) компоненты поддерживающей силы согласно 
(10), будем иметь 


„о Пе 40, О=Олб[у (946 (15) 

















Ю = Р.В", 5т(®- 1)2п $т(® + 12 2, 
° 4В ®-—1 ©-+1 ' ор. 
| И то, — } 6. 
| 2В 2 ®-—1 ©+| 5— 
4— 
* 2 1 1 3 
В = р ‚ В= № 8. и 
28-" 26 я 
(16) :] 
Численный анализ аналитических выраже- 


ний (15) и (16), приведенных на рис.2, показывает: 
1. При с«=1/2 условие периодичности гидро- в я: а | 

динамического давления (условие замкнутости сма- НЫ 

зочного слоя) выполняется. При этом значении ®« — рРис.2. Зависимость безразмерной несущей способ- 

несущая способность подшипника на 40-50% выше, = ности от параметров | и сэ при различных значени- 

чем при 0. ях параметра связи № и параметра сжимаемости Л: 
2. С увеличением значения параметра связи 1 - №=0,8; Л=100; 2 - М№=0,8; Л=200; 

М и параметра сжимаемости Л несущая способность — 3-№=0,95; Л->оо 

подшипника возрастает. Наиболее оптимальными по 
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несущей способности и силе трения значения параметра № принадлежат интервалу 
№ Е[О,95; 0,98]. 

3. При ®«=1/2 предложенная конструкция подшипника с адаптированным профилем опор- 
ной поверхности имеет не только повышенную несущую способность, но и меньшие габариты, по 
сравнению с подшипником с круговой (‹«=0) опорной поверхностью. 

Выводы. Предложено точное автомодельное решение задачи о гидродинамическом расчете ра- 
диального подшипника, работающего на сжимаемой микрополярной смазке. Найдены оптималь- 
ные по несущей способности подшипника значения параметров, присущих микрополярным смаз- 
кам. Теоретически обоснован профиль подшипника, обеспечивающий повышенную его несущую 
способность. 
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